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We use the optimized perturbation theory, or linear 5 expansion, to evaluate the critical exponents 
in the critical 3d O(N) invariant scalar field model. Regarding the implementation procedure, this 
is the first successful attempt to use the method in this type of evaluation. We present and discuss 
all the associated subtleties producing a prescription which can, in principle, be extended to higher 
orders in a consistent way. Numerically, our approach, taken at the lowest nontrivial order (second 
order) in the S expansion produces a modest improvement in comparison to mean field values for 

Othe anomalous dimension n and correlation length v critical exponents. However, it nevertheless 
points to the right direction of the values obtained with other methods, like the e-expansion. We 
. discuss the possibilities of improving over our lowest order results and on the convergence to the 

known values when extending the method to higher orders. 
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I. INTRODUCTION 



-o - 

A distinct feature as we approach second order phase transition points is the emergence of critical phenomena 
and associated universality and scaling properties, as a result of diverging correlation lengths. Thanks to these 
O , characteristics we can relate many different systems just by simple general distinctions, like dimension and symmetry, 
with their universal behavior set by critical exponents, independent of the microscopic dynamics 0. On the other 
| hand, close to the critical temperature of transition we are faced with the problem of breakdown of perturbation 
^ i methods in field theory due to the appearance of infrared divergences, as the correlation lengths diverge (or masses 
' vanish), which then require the use of nonperturbative methods to study the physical system around the critical points 
and, ultimately, the phase transition process itself. The present work was motivated by the recent progress in dealing 
with field theory phase transitions in the context of the nonperturbative method of the optimized mass, or linear 5 
expansion (LDE) HHSI3. 

Here, for the first time, this method is used in the derivation of the critical exponents 
■ of the three-dimensional O(N) invariant A</> 4 model. 

Critical exponents have been extensively studied at different levels of accuracy, particularly within the e-expansion 
and renormalization group (RG) methods (see, e.g. Ref. for a throughout discussion of various methods of 
evaluating critical exponents in field theory). However, it is well known that these methods may not be applicable 
or appropriate for a number of important physical systems. Examples include, for instance, the reliability of the 
e-expansion in the study of the electroweak phase transition, where there may not be a visible fixed point in three 
dimensions when we perform an expansion around dimension d = 4 — e 7]. There are also the important cases of 
multicritical phenomena in two dimensions with no obvious upper critical dimension or with no O(N) invariants, in 
which case nonperturbative techniques like the e-expansion and large- N approximation break down The use of 
renormalization group techniques have shown recently to produce impressive good results to critical exponents [jj, 
however, there we commonly have to deal with a set of flow equations that may become quickly very complicate to 
solve for complex systems. 

Using the LDE one makes use of an interpolation of the original action by modifying the masses and coupling 
constants in terms of a fictitious parameter 5 that at the end is taken to the unity, so as to recover the original model. 
At the same time an arbitrary mass parameter must be introduced to balance the dimensions of the interpolated 
theory. At the end it is fixed in such a way that nonperturbative information can be taken into account in quantities 
which have been perturbatively computed. Concerning the O(N) invariant scalar field model in three dimensions, the 
convergence properties of the LDE have been studied recently in great detail, specifically at the critical point of phase 
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transition in the large- N limit |j, |g, [l(J . Those results are also encouraging concerning the computation of critical 
exponents as performed in this paper. An advantage of the LDE method when compared to other nonperturbative 
methods is that we remain all the time within the familiar grounds of perturbation theory. The method then becomes 
of much simpler use than the traditional ones and we expect that the LDE may even further improve the accuracy 
of other nonperturbative methods when used in conjunction with them. We also expect the method to be suitable 
in the evaluation of the critical properties in the transition point in those cases where no clear expansion parameter 
exists or when the conventional techniques do not apply, as in the models discussed in the previous paragraph. The 
LDE has been successfully applied before to study field theory at finite temperature and phase transitions in those 
situations where conventional perturbation theory breaks down 0,0, El- I n particular, in Ref. 12], the LDE was 
applied to study the phase transition patterns in coupled scalar field models reproducing previous renormalization 
group applications to the problem and also predicting new critical points, showing the advantages of the method for 
applications regarding complex systems that may exhibit multicritical points. 

In this work we shall consider the same O(N) model studied in Refs. 0, H S Q. It is described by the 0(N) 
invariant super-renormalizable action 



(1.1) 



where ro = r + A, with A representing the mass renormalization counterterm needed to remove any ultraviolet 
divergence. The full propagator for this theory is given by 



G(p) = [p 2 + m 2 + E rcn (p) - E ren (0)] 1 , (1.2) 

where m 2 = r + E ren (0) with £ rcn representing the renormalized self-energies. At the critical point, m 2 — 0, which 
implies the Hugenholtz-Pines theorem, r c = — S rcn (0). Using the model described by Eq. (|1 . 1|) is advantageous since 
our results can readily be compared with other results that have been extensively obtained in connection with the 4 
O(N) invariant scalar model. This is a particularly important model also for the reason that it can describe the critical 
properties (due to universality) of many different physical systems, like polymers (N = 0), Ising models (N = 1), 
superfluid Helium, Bose-Einstein condensation of atomic atoms (TV = 2), ferromagnets (N = 3), the scalar sector of 
the electroweak standard model (N = 4), etc. A similar attempt to compute critical exponents with the S expansion 
was performed by the authors of Ref. |13j . However, they made use of the nonlinear version of method which leads 
to a considerable more complicated perturbation series (see also 0]), which quickly becomes cumbersome beyond 
leading order in S and their results were not so good for the same critical exponents evaluated here. On the other hand, 
the way the linear version is employed here avoids those difficulties and the results can in principle be extended to 
arbitrary orders as usual perturbative calculations. Our results for the correlation length (or mass) critical exponent v 
and the anomalous dimension exponent 77, already to the lowest nontrivial order (second order) in S, though showing 
only modest improvements over the mean field values, are shifted away from these in the direction of the values known 
from high precision numerical results despite the simplicity of our calculations. 

This paper is organized as follows. In Sec. II we briefly review the LDE method and present the interpolated version 
of the O(N) invariant scalar field action relevant to our study. In Sec. Ill, we carry out the formal evaluation of the 
critical exponents for the correlation length v and anomalous dimension r\. Our concluding remarks arc presented in 
Sec. IV, where we also discuss the extension of our method to higher orders. 



II. LDE AND THE INTERPOLATED EFFECTIVE SCALAR MODEL 



Let us start our work by reviewing the application of the LDE method to our problem. The LDE was conceived 
to treat nonperturbative physics while staying within the familiar calculation framework provided by perturbation 
theory. In practice, this can be achieved as follows. Starting from an action S one performs the following interpolation 



S -> Ss = SS + (1 - 6)S (M) , (2.1) 

where So is the soluble quadratic action, added by an (optimizable) mass term M, and 5 is an arbitrary parameter. 
The above modification of the original action somewhat reminds the usual trick consisting of adding and subtracting 
a mass term to the original action. One can readily see that at 6 = 1 the original theory is retrieved, so that S actually 
works just as a bookkeeping parameter. The important modification is encoded in the field dependent quadratic term 
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Sq(M) that, for dimensional reasons, must include terms with mass dimensions (M). In principle, one is free to choose 
these mass terms and within the Hartree approximation they are replaced by a direct (or tadpole) type of self-energy 
before one performs any calculation. In the LDE they are taken as being completely arbitrary mass parameters, which 
will be fixed at the very end of a particular evaluation by an optimization method. One then formally pretends that S 
labels interactions so that Sq is absorbed in the propagator whereas SSq is regarded as a quadratic interaction. So, one 
sees that the physical essence of the method is the traditional dressing of the propagator to be used in the evaluation 
of physical quantities very much like in the Hartree case. What is different between the two methods is that, within 
the LDE the propagator is completely arbitrary, whereas it is constrained to cope only with the so-called direct terms 
(i.e. tadpoles) within the Hartree approximation. So, within the latter approximation the relevant contributions are 
selected according to their topology from the start. 

Within the LDE one calculates in powers of S as if it was a small parameter. In this respect the LDE resembles 
the large- N calculation since both methods use a bookkeeping parameter which is not a physical parameter like the 
original coupling constants and within each method one performs the calculations formally working as if N — ► oo 
or 6^0, respectively. Finally, in both cases the bookkeeping parameters are set to their original values at the end 
which, in our case, means 5 = 1. However, quantities evaluated at any finite LDE order from the dressed propagator 
will depend explicitly on M, unless one could perform a calculation to all orders. Up to this stage the results remain 
strictly perturbative and very similar to the ones which would be obtained via a true perturbative calculation. It is 
now that the freedom in fixing M generates nonperturbative results. Since M does not belong to the original theory 
one may require that a physical quantity calculated perturbatively to order-<5 fc be evaluated at the point where 
it is less sensitive to this parameter. This criterion, known as the Principle of Minimal Sensitivity (PMS), translates 
into the variational relation Hal 



dM 



. 



M,S=l 



(2.2) 



The optimum value M which satisfies Eq. 11'. 21 must be a function of the original parameters including the 
couplings, which generates the nonperturbative results. 

Following the considerations above, in order to interpolate the original theory described by Eq. one may use 



So = x[|V0| 2 +r4> 2 +M 



as in Refs. 3]. Then, the interpolated action reads 



2 J.21 



(2.3) 



dfx 



I| V 0| 2 + i(r + M 2 ) 



-M 



Su 

IT 1 



-M4> 2 



(2.4) 



where As represents the renormalization mass counterterm for the interpolated theory, which depends on the param- 
eters M and 5. It is important to note that by introducing only extra mass terms in the original theory the LDE 
does not alter the polynomial structure and, hence, the renormalizability of a quantum field theory. In practice, the 
original counterterms change in an almost trivial way so as to absorb the new M and 8 dependence (for details, see 
for instance 0). The compatibility of the LDE with the renormalization program has been shown in the framework 
of the O(N) scalar field theory at finite temperatures, in Ref. showing that it consistently takes into account 

anomalous dimensions in the critical regime. 

For the interpolated theory, the full propagator G^ (p) , can be written as 



(2.5) 

where Sren are the (renormalized) momentum dependent self-energies evaluated in powers of 6 and to 2 = r + Erfn (0) . 
At the same time, the bare propagator is given by 



G^(p) = p 2 + m 2 + (1 - 5)M 2 + EW (p) _ E W (o) 



G^(p) = \p 2 + n 2 }- 1 , (2.6) 

where il 2 = rn 2 + M 2 . It is interesting to note that at the critical point (to = 0) the propagator given by Eq. I|2.(j|) 
does not generate any infrared divergences since it is automatically regulated by the LDE mass parameter M whereas 
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the equivalent bare propagator of the original theory is massless. In fact, this is the main problem concerning any 
eventual perturbative evaluation of quantities at the critical point, an example being the problem of calculations of 
critical temperature shifts, AT C 0,U0,0|- 

In our evaluation of the renormalized quantities and fixed points we shall use the following Feynman rules for the 
vertices: <5M 2 , — SA$ for the quadratic interaction and —Su for the quartic one. Those rules are used to evaluate, 
perturbatively in 5, all the relevant physical quantities. After that, as discussed in Sec. II, the nonperturbative results 
are produced by means of the PMS condition, Eq. (|2.2|) . 



III. CRITICAL EXPONENTS EVALUATION UP TO ORDER-d 2 



We now proceed to the evaluation of the critical exponents for the correlation length and anomalous dimension, 
which follow from the usual scaling relations at the critical point. Following the standard conventions and definitions 
given, e.g., in Parisi's textbook one writes the scaling relations as 



oc (r - 5r s c f Z 2 = 



dr 
dm 2 



Z-i m 



2C 2 



C2 = — - 1 
2v 



(3.1) 



and 



\ dG^\p) 

{ dp* 



P 2 =o 2 



(3.2) 



where Zi is the field renormalization function and A is the effective renormalized coupling, defined in terms of the 
one-particle irreducible four-point function T^O) as 



\ = Zf 



2 r 4 (o) 



(3.3) 



In order to compute the critical exponents defined in Eqs. (|3.1|l and (|3.2() . we need to compute the (3- function, 



'-■"bp* --§£*• (3 - 4) 

from which the fixed points are obtained. Note that g — u/m is taken as a dimensionless coupling. The fixed points 
are then obtained, as usual, from the solutions of (3{\ c ) — 0. 
From the relations (|3.1|) and l|3.2|) we can also define 



CMA^-f-ln^A)], (3.5) 

and 



C 2 (A)=-|— ln^ 2 (A)/^(A)], (3.6) 

which, when evaluated at the fixed point, will determine the constants C\ = Ci(A c ) e C2 = C2(A C ) defining the critical 
exponents in l|3.1|l and l|3.2|l . 

Having introduced the quantities given above, let us now use the LDE to evaluate r^O), Zi and the self-energy 

Srfn expanded in powers of <5. Note that at first order in 8 the results are trivial, since it corresponds just to the 
lowest order one-loop perturbative expansion for which Z\ = Z% = 1. In this case, using the definitions (|3.5|l and 
(|3.6() . one obtains C\ = C2 = which lead to the well known mean- field critical exponents 77 — and v = 1/2. 
The next order results in 6 will however already lead to results departing from the mean-field ones. For consistency 
we carry the calculations of all quantities to the first non-trivial order, 5 2 . Also, all our calculations are performed 
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using dimensional regularization and the modified minimal subtraction (MS) renormalization scheme which amounts 
to replace the the momentum integrals by 



f d 3 p f _ ( £*jg Y f d 2uj p 
J (2tt)3 ~* J p - { 4vr ) J (2tt) 2 - ' 

where 2lo — 3 — 2e while u s is an arbitrary mass scale and je — 0.5772 is the Euler-Mascheroni constant. 



A. Evaluating the optimized self-energy. 



As already emphasized, all relevant physical quantities evaluated perturbatively with the LDE will depend on the 
arbitrary quantity Q 2 = m 2 + M 2 , which is present in the LDE bare propagator, Eq. I|2.6|l . Nonperturbative results 
can be generated by optimizing the renormalized self-energy, Sren(O), which then generates the optimum M. Having 
mass dimensions, this quantity can be expressed as a function of the original m and u parameters which, as we shall 
see, will allow us to express the optimum quantity f2 in terms of m and the dimensionless coupling g = u/m. The 

choice of Eren(0) as the physical quantity to be optimized is well justified since this is an important quantity in the 
study of critical phenomena signaling phase transitions via the Hugenholtz-Pines theorem. It also fixes the effective 
scale determining all critical exponents through the scaling relations at the critical point. To order-<5 2 the self-energy 
has been explicitly evaluated in Ref. When m ^ it is given by 1 



1(2) 



(0) = - 



Suil { N + 2 



8tt V 3 
5 2 u 2 (N + 2) 



6 2 uM 2 fN + 2 



(8tt) 2 18 

We then set S = 1 solving the PMS equation 



16ttS1 

MS- 



3 

2 + 41n 



2„,2 



S u 

1287T 2 



N ■ 



0(S 3 ) 



(3.8) 



dM 



_=0, 

M 



(3.9) 



to obtain the (iV-independent) roots 



Mo = 
± (12tt) : 



{u 2 ±u[(l2irm) 2 + u 2 ] 1/2 } 



(3.10) 



The solutions Mq and M_, at the critical point m = 0, are trivial ones, while M + (m —►{)) = u/(6ir) remains nonzero 
even at the critical point and it is then able to effectively lead to nonperturbative results as shown below, apart from 
agreeing with our previous results Q . This nontrivial root (omitting the index "+" from now on) can be written in 
terms of the dimensionless coupling, g = u/m, as 



—2 m 

M = 



72tt 2 



g 2 + Ung 



-,1/2- 



(127T) 



(3.11) 



2 2 

which is, manifestly, a nonperturbative quantity. Finally, the optimized Q — m 2 + M can be written as 



= m 2 \l + ^[g 2 + 12ngF(g)]\ , (3.12) 



1 See Ref. for details concerning the m = case. 
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where we have defined 



1/2 



F(g) = 



(12vr) 2 



(3.13) 



So, in the following, our strategy will be to evaluate all relevant quantities with the LDE to order-5 2 and then perform 
the substitution fl — > f2 where the latter quantity is given by Eq. I|3.12[l . 



B. Evaluating the critical coupling g c 

Let us start by evaluating the (3 function whose roots define g c . The relevant Green's function is the four-point one, 
which to order-<5 2 reads 

3 r2 2 {N + 8\ f d 3 q 1 



r 4 (0) = Su - -S 2 u 



2 V 9 ) J (2tt) 3 (q 2 +ti 2 ) 2 
= Su-^(^), (3.14) 

where the first term on the RHS is the tree vertex and the second term is the one-loop correction. In addition to 
the four-point function, according to Eq. QH.3JI . one needs the field renormalization function in order to define the 
effective renormalized coupling. To evaluate Zi, as given by Eq. 13. 2|) . we need the full propagator G^(j>) which, to 
order- J 2 , is 

G^(p) = [p 2 + m 2 + Eg) (p) S^(O)]- 1 , (3.15) 
where E^(p) - £&(0) is given by E3 



vWrs v(2), n s (N + 2)S 2 u 2 \^ 3fl fp\ l/p 2 + 9n 2 \ 

E&(p) " (0) = - 18(47r)2 [l - - -tan (_) - - In j 

Substituting Eqs. (|3.16[) and l|3.15|l in (|3.2I) . one gets 



(3.16) 



which, after expanding to order-5 2 , becomes 



18(47r) 2 540 2 



Zl 1 ~ 18(47r)W ■ (3 - 18) 



From Eqs. (|3.3|l . (|3.14|l and Ij3.18|) we then obtain the renormalized coupling to order S 



2 



*=--^(^)+0(P) , (3.19) 



77) 



IQirmQ, \ 9 



from which the 0- function will follow via Eq. (|3.4|) . But before performing the derivatives with respect to g (or m) 
let us recall that this quantity also appears in f2. Then, substituting the optimized quantity Eq. I|3.12|l in H3.19fl and 
setting 5 — 1 we obtain the optimized renormalized coupling given by 

A(g) = g g2( f g + 8) {l + ^ [g 2 + 12ngF(g)]y 1/2 , (3.20) 
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which is, of course, a nonperturbative quantity. Now, by using the definition for the /3-function, which from Eq. (|3.20() 
can be expressed as 

<•■»> 

one finds the fixed points, g c , as given by the solutions of (3(g c ) = 0. For a fixed value of N ^ we always find three 
roots, a trivial one (g c = 0), a purely imaginary one (which leads to unacceptable scaling relations at the critical 
point) and a positive real root. This last one can easily be found numerically, and is given, e.g. for a few cases, by 
g c {N = 1) = 15.3524, g c (N = 2) = 12.2968, g c {N = 3) = 10.3692. The nontrivial real positive roots are the ones that 
will be used in the subsequent determination of the critical exponents. 

C. Evaluating the critical exponents 

Now, the evaluation of the critical exponent r\ is fairly easy. From Eqs. I|3.2[) and (|3.5|l we have that 

f] = -g^-\n[Z 1 (g)} , (3.22) 

Og g=g c 

where Z\ represents the optimized field renormalization function. This quantity can be obtained directly from Eq. 
(|3.18|l with the replacement — ► fl. Then, setting 5=1, one obtains 

^= i -f»{ i+ 7^[» 2+i2 ^r 1 ' < 3 - 23 ' 

From Eq. 13.22fl and the previous results for the fixed points, we then obtain the results (for some representative 
values of N) tj(N = 1) = 0.0026, r](N = 2) = 0.0029, r](N = 3) = 0.0030. These values are contrasted to the results 
from other methods in Table I. Results for other values of N can also be easily obtained from the previous equations. 

We now turn to the calculation of the critical exponent v to order-^ 2 . According to our prescription, to evaluate 
this quantity one computes the optimized Zi which is given by 

Z 2 = -^, (3.24) 

(2) (2) 

where r = m 2 — S ron (0). The optimized self-energy at order 5 2 , E rcn (0), can be trivially obtained from Eq. i|3.8|) by 
performing the substitutions M — > M and £1 — > fl as given by Eqs. (|3.11(l and 13.12J1 . The critical exponent v then 
follows from Eqs. (|3.6|) and (|3.1|) . or 

- g—HZ2(9)/Zi(g)}\ , (3.25) 

og ) g=ga 

which, at 8 = 1 and for a few representative values of N, yields the results v(N = 1) = 0.5287, v(N = 2) = 0.5362, 
v(N = 3) = 0.5422. These results, together with those for rj, are shown in Table I so that they can be compared with 
the results obtained using different approximations. 

The e-expansion results for the critical exponents can be found e.g. in Ref. 0. In Parisi's book ^(| one can find, for 
instance the results for the critical exponents for one and two-loop perturbation theory (PT). In Ref. it is used 
a variational perturbation theory (VPT) to continue the renormalization constants of three-dimensional </) 4 theories 
to the regime of strong bare coupling. The authors of Ref. 0] have used the logarithmic (^-expansion (DE-log) to 
obtain the exponents rj and v. All the results from these different methods are shown in Table I together with ours. 

IV. CONCLUSIONS 

In this paper, the LDE has been applied for the first time in the evaluation of critical exponents. An early attempt 
to use this type of method was carried out by Gandhi and McKane 01 wn0 used a variant of the method, known as the 
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TABLE I: Numerical results for v and 77. 



Method 


(N=l) 
v (N=2) 
(N=3) 


(N=l) 
r, (N=2) 
(N=3) 


LDE - 0(5 2 ) 


0.5287 
0.5362 
0.5422 


0.0026 
0.0029 
0.0030 


e — expansion 


0.626 
0.655 
0.679 


0.037 
0.039 
0.039 


PT one — loop 


0.6 
0.625 
0.647 







PT two — loops 


0.630 ± 0.015 


0.031 ± 0.004 


VPT 


0.63 
0.670 
0.706 


0.030 
0.032 
0.032 


DE - log 


0.53 


0.0013 



logarithmic 5 expansion. However, the authors did not implement the method correctly, since in their implementation 
they have not made use of any optimization scheme, just plain expansion in <5, which led to many difficulties associated 
to the fact that their fixed point coupling g c behaves like 1/5. 

Here, our main concern was to show how to circumvent the difficulties found in Ref. 01 . Our approach to the 
problem is based on the central idea of first obtaining optimized values for the arbitrary parameter M by applying the 
optimization procedure PMS to the renormalized self-energy at a given order in S, Sffn(0), which then generates an 
optimum value for M in terms of the original parameters of the theory (m and g — u/m). Then, the full propagator 
and four-point functions are evaluated to the same perturbative order with the LDE and their optimum values obtained 
by performing the substitution M — > M(m,g). Our analytical results show explicitly that the M dependence on g is 
nonperturbative (all orders are present). The relevant optimized quantities Z\, Zi and the (3 function are obtained 
upon using derivatives with respect to g. This whole procedure then generates the optimized values for the fixed point, 
<? c , as well as for the critical exponents r\ and v (note that, in contrast to the approach used in |13|. the optimization 
procedure is already performed prior to computing the /^-function and the fixed points). In principle this prescription 
can be extended to any perturbative order in 5. The only difficulties of extending our results to higher orders are 
technical ones, associated to the evaluation of multi-loop diagrams containing massive propagators. 

As far as our numerical values are concerned, they already at nontrivial lowest order show improvement, though 
modest, over the mean field theory ones. In a sense, the physical quantities evaluated by us with the LDE to order-5 2 
contain only two-loop diagrams so it comes as no surprise the fact that our lowest order results are not as good as 
the ones obtained by resumming higher order corrections. One expects, as shown by the many LDE applications, 
that by considering higher order contributions one can quickly obtain better results from a calculation which has the 
advantage of being completely perturbative regarding the evaluation of Feynman diagrams. In this context, from our 
previous experience 01- Q in extending to higher order the LDE to obtain the a the critical temperature shift for 
an homogeneous Bose gas, which is a problem that share many similarities to the one studied here, we can advance 
a number of important issues that must be handled. First is the slow convergence behavior observed in 0] by 
the direct extension of the LDE method studied here. Since both the calculation of the temperature shift for Bose 
condensation (for N = 2) and the exponents computed here make use of the same critical model, it is reasonable to 
expect that the same convergence behavior will show up in higher order calculations of the critical exponents within 
the LDE method. The general explanation and throughout understanding of the slow convergence behavior for the 
Bose condensate critical temperature shift was given in Ref. p|. Also, in the same reference, it was discussed how 
appropriate resummation techniques can speed convergence. The same resummation techniques used there could also 
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be used here in order to improve convergence. A second important issue that commonly arises when extending the 
LDE beyond second order is the appearance of imaginary solutions upon the use of the optimization, Eq. 12. 2|) . In 
previous references we have shown that only the real part of those solutions make physical sense and indeed they lead to 
correct results, despite some obvious embarrassment of having to deal with those imaginary solutions. To circumvent 
this problem an alternative optimization has been proposed recently |l8j . which make use of additional parameters 
from third order onwards. In practice this can be achieved by replacing, in Eq. (|2.1|l . terms like M* 2 = (1 — S)M 2 by 
terms like 



M* 
M* 



[l - a x 8 - (1 - ai)<5 2 ] 1/2 M, at third order 



1 



ai5 - (1 - ai)<$ 2 + ajS 2 {l - 6)] 1/2 M, at fourth order 



(4.1) 



and so on, where the a,-, j = 1, . . . , n — 2 are additional parameters added at order n in the LDE and to be determined 
with M by generalizing the optimization condition Eq. (|2.2|l to a system of equations, 



dM 



,8=1 



= o, 



dM 2 



M,a\,.. .,a n — 2 ,5=1 



o, 



jra-l<jj(n) 

=o (4 2) 

dM 71 ' 1 M,a,...,a„-2,S=1 ' 

As shown in |18| this optimization procedure generates, in the Bose-Einstein condensation case, only real values for 
the LDE interpolating parameters with the added bonus of fast convergence already at lowest orders. As far the LDE 
application of Ref. [18J to the Bose-Einstein condensation problem is concerned, it has produced some of the most 
precise and stable analytical predictions for the critical temperature shift. Those results thus also give an indication of 
the applicability of the use of the LDE in analogous models close and at the critical point. We are currently working 
on the use of this modified optimization to the problem studied in this paper and we will report the results elsewhere. 
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